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ABSTRACT 

The line of sight to the quasar B0218+357, one of the most studied lensed systems, intercepts a z = 0.68 spiral galaxy, 
which splits its image into two main components A and B, separated by ca. 0.3", and gives rise to molecular absorption. 
Although the main absorption component has been shown to arise in front of image A, it is not established whether some 
absorption from other velocity components is also occuring in front of image B. To tackle this question, we have observed 
the HCO + (2-l) absorption line during the commissioning phase of the new very extended configuration of the Plateau 
de Bure Interferometer, in order to trace the position of the absorption as a function of frequency. Visibility fitting of 
the self-calibrated data allowed us to achieve position accuracy between ~ 12 and 80 mas per velocity component. Our 
results clearly demonstrate that all the different velocity components of the HCO + (2-l) absorption arise in front of the 
south-west image A of the quasar. We estimate a flux ratio Ja/Jb = 4.2^ 1 'o at 106 GHz. 
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1. Introduction 

Molecular absorption at intermed iate redshifts has been de- 
tected only in a f ew objects C e.g. I Wiklind fc Combesl llflM 
119961 H997Lll998tlKanekar et al.ll2005lT Among these, three 
absorption systems are caused by a galaxy lying on the 
line-of-sight of a quasar and acting as a gravitational lens. 
Two systems, one in front of PKS1830-211 (at z = 0.89), 
the other in front of B0218+357 (at z — 0.68) are detected 
at millimeter wavelengths in the lines of several molecules, 
including HCO+, HCN and H 2 CO. 

The light rays associated with the quasar PKS1830— 211 
form two bright images that probe different regions of 
the intervening galaxy. This provides information on the 
latter's kinematics, its mass, as well as on the physi- 
cal and chemical conditions in its interste llar medium 
(| Wiklind fe Combed H99l iMuller et aHliMOll . These two 
papers showed that the two z ~ 0.89 absorption compo- 
nents of PKS1830— 211 are associated each with one of the 
two gravitational images. 

The second mm-absorption system occurs in the line 
of sight to B0218+357 and is caused by a galaxy at a 
redshift of z = 0.68466 (TBrowne et all 119931 : ICarilii et all 
Il993t IWiklind fc Combes! 1 19951) . Two bright images of the 
quasar, hereafter referred to as A (to the SW) and B 
(to the NE), have been resolved at rad io cm wavelengths 
jO'Dea et al.lll992HPatnaik et al.lll993f) . The distance AB 
between the two images is ~ 0.3", the smallest angular 
separation among the known galaxy-mass lenses. The flux 
ratio is Ja/Jb ~ 3 between 5 GHz and 22 GHz. Image B 
lies in the center of an Einstein ring whose diameter is sim- 
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ilar to the distance AB. Each image reveals intricate sub- 
struct ures at very high angular resolution (|Patnaik et al.l 
119951: iBiggs et al]|2003t ). The constraints provided by the 
complex image pattern and by the time variability of 
the background source flux, make of B0218+357 one of 
the best objects to measure Hp at intermediate redshift s 
(IBiggs et alil2003t IWucknitz et alJl2004t lYork et al"1l2005lb 
Deep ACS/HST observations of B0218+357 bv lYork et al.l 
(2005) reveal that the lensing object is a spiral galaxy seen 
nearly face-on and whose center lies close to image B, the 
center of the Einstein ring. Accordi ng to the lensing model 
proposed bv lWucknitz et all (|2004l ) . image A is located at 
about 2 kpc from the center of the lensing mass distribu- 
tion, and image B at 0.4 kpc. HI absorpt ion is detected over 
a velocity width of a bout 100 km s _1 (ICarilii et alj 1993; 
i Kanekar et al.l 120031 ). in front of image A ( Carilli et al.l 
although so me absorption may a lso occurs in front of 
the Einstein ring (jKanekar et al II2003T). H 2 CO (2 12 - 2 n ) 
absorption has been observed bv lMenten fc RekH (|1996| ) to- 
wards image A, with a total width of ~ 12 km s _1 , much 
narrower than the HI profile. Because of limited signal-to- 
noise ratio, they could however not exclude the possibility 
of absorption in front of B. 

FigQ] shows the HCO + (l-2) absorption profile towards 
B0218+357 observed with a high sensitivity and velocity 
resolution (Muller et al., in prep.). Like the HCN (1-2) pro- 
file, it shows at least four velocity components over a width 
of 25 km s _1 , much broader than that of H 2 CO (2 i2 — 2n). 
The component with the deepest absorption matches in 
width and velocity the H 2 CO (2i 2 — 2n) absorption profile 
and most probably originates like this latter, in front of im- 
age A. The other three components, however, had not been 
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observed so far with a high angular resolution and their 
location, in front of A, B, or of the Einstein ring remained 
unknown. Yet, the knowledge of this location is essential 
for the interpretation of the molecular absorption data in 
terms of line opacities, molecular column densities and gas 
kinematics. 

In this letter, we present observations of the HCO + (l-2) 
line toward B0218+357, obtained with the new, extended 
configuration of the Plateau de Bure Interferometer, which 
allow us to trace the position of all four absorption compo- 
nents. 

2. Observations and data analysis 

In the frame of a survey of molecular absorption lines to- 
wards B0218+357, we have observed the HCO+(l-2) line 
in June 2005 and July 2006, with a compact configura- 
tion of the IRAM Plateau de Bure Interferometer (PdBI). 
These observations were self-calibrated on the continuum 
source, which was not resolved with the ~ 4" FWHP syn- 
thetized beam. The resulting average HCO + (l-2) spectrum 
is shown in FigQ] The absorption extends over ~ 25 km s _1 
and can be decomposed into four Gaussian velocity com- 
ponents, with FWHM of 4 - 5 km s -1 . 

We have further observed on November 8 th 2005 the 
HCO + (l-2) line with an extended configuration (maximum 
projected baseline B max of 385 m) and on November 9 th 
2005 with the new, very extended configuration of the PdBI 
(B max =650 m). We emphasize that these observations were 
the first science data acquired with the new extended con- 
figuration and that they were performed during its com- 
missioning phase. In particular, the instrument baselines 
were not yet properly calibrated, the on-source integration 
time was limited (110 min in each configuration), and the 
weather not at its best (rms phase deviation prior to self- 
calibration of ~ 50° on the longest baselines). Therefore, 
the data needed to be self-calibrated on the continuum 
source, which was half-resolved in the EW and NEW di- 
rections. 

As in the June and July sessions, we observed simultane- 
ously a 540 MHz-wide frequency band centered at 105.690 
GHz, consisting of two overlapping 320 MHz-wide con- 
tinuum sub-bands (L03 and L04) with a channel spacing 
of 2.5 MHz (7.1 km s" 1 ), and a 80 MHz-wide sub-band 
(L01) centered on the redshifted HCO + (T2) line frequency 
^=105.882 GHz, with a channel spacing of 0.31 MHz (0.9 
km s _1 ). The high frequency edge of the continuum band 
was dropped and its bandwidth restricted to 440 MHz in 
order to avoid contamination by the HCO + (l-2) line. 

The method used to self-calibrate on the half- 
res olved continuum sou rce was similar to the one used 
by iMuller et al.l (|2006[ ). for localizing the absorption 
components in front of PKS1830— 211. The data from 
November 8 th and 9 th were calibrated separately, using the 
GILDAS/CLIC software, with the following procedure: In 
a first step, the radio frequency (RF) bandpass of L03 and 
L04 was calibrated directly on the B0218+357 continuum 
signal. The main purpose of using these absorption-free 
bands was to correct for short term phase and intensity vari- 
ations linked to atmosphere fluctuations by self-calibrating 
on the continuum. L01 was calibrated in RF separately, us- 
ing low order polynomials (first degree in amplitude and 
third degree in phase). The narrowness of the absorption 
line (< 10 MHz) with respect to the L01 bandwidth (80 



MHz) ensured that the calibration was not badly affected 
by the absorption features. Note that we choose a baseline- 
based RF calibration, rather than an antenna-based, as the 
signal intensities for the longest baselines were severely af- 
fected by atmospheric phase decorrelation. 
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Fig. 1. Spectrum of the HCO + (l-2) absorption towards 
B0218+357 as observed with the PdBI in a compact config- 
uration. Gaussian fits of the different velocity components 
a, b, c, d are overlaid as dashed lines. 

Next, the continuum sub-bands were self-calibrated in 
amplitude and phase by calculating the complex gains cor- 
responding to a system of two point-like sources separated 
by AR.A. = 309 mas, ADec. = 128 mas, and with a flux 
ratio of 4.2 (a justification of theses values will be given in 
§3, and/or can be found in the reviews by IWucknitz et alj 
120041 : iMittal et al.ll2006a|bl ). The total continuum intensity 
(~ 0.4 Jy) was normalized to unity. The gains calculated for 
the continuum sub-bands were then applied to the visibil- 
ities of the L01 sub-band channels. Finally, the continuum 
visibilities, as derived from the double-source model, were 
subtracted from the calibrated L01 visibilities. 

The calibrated L01 datasets from November 8 th and 
9 th were then combined. The visibilities corresponding to 
uv radii larger than 200 m were fitted, channel by channel, 
with a single point-like source representing the position and 
strength of the absorption signal. Visibilities with shorter 
uv radii were discarded, as they bring little information on 
the signal position. The fitted source intensities showed a 
small offset (~ 10% of the total continuum) with respect 
to zero (the continuum level as defined from the continuum 
sub-bands) - probably the result of residual IF and/or RF 
bandpass calibration errors. This was taken care of by av- 
eraging the visibilities of the L01 channels that are free of 
any absorption and by subtracting the so-computed resid- 
ual continuum visibility from all L01 channels. The fit of 
the point-like absorption sources was then repeated. 

Fig[2]shows the result of the final fits. The upper plot (i) 
shows the absorption profile as a function of velocity, while 
the two middle plots (ii and Hi) show the position of the 
absorption (or, more exactly, the position of the centroid 
of the absorption) as a function of velocity. The removal of 
the continuum insures that zero intensity (0) and zero offset 
(0,0) effectively correspond to the intensity and position of 
the centroid of the continuum sources. The source position 
accuracy in the individual velocity channels was calculated 
by the GILDAS fitting routine UVFIT. It can be crudely 
expressed as ~ 6>b ca m / 2 SNR, where #bcam is the FWHP 
of the synthesized beam (~ 1.2 x 0.7", P.A. 176°) and SNR, 
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the signal-to-noise ratio on the source intensity. The noise 
level, per 0.9 km s _1 velocity channel, is 3.4% of the con- 
tinuum intensity. For all channels with absorption signal 
(—10 < Vhel < +12 km s _1 ), the fitted phases indicates 
that, within the uncertainties, the absorption arises from a 
small size region that coincides in position with image A. 
We further averaged the position offsets per velocity com- 
ponents, with weights equal to the square of the intensity 
of the Gaussian profiles. The average positions are given in 
Table H and Fig©?;. 

3. Discussion 

Our data show that all the molecular absorption com- 
ponents (from Vhel = —10 to +12 km s" 1 ) originate 
in front of image A. This is not surprising since im- 
age A appears more obscu r ed in the visible than im- 
age B ( Grundah l fc Hiorthl [l995). As was mentioned, 
IMenten fc Reidl (|l996l ) have already shown that the Vhel 
~ 5 km s _1 H2CO absorption component arose in front of 
image A. 

The background radio source presents a core-jet mor- 
phology that is reflected in both lensed images in the form 
of sub-structures. Those consist mainly in two main com- 
ponents, labelled Al and A2 (and Bl and B2 ), which are 
separated by ~ 1.4 mas (Pat naik et al.|[T995f h High sen- 
sitivity radio imaging reveals in addition the presence of 
a weak, relatively exte nded (10-15 mas) knotty jet and of 
a possible counterjet (|Biggs et al.l [2003). Our source po- 
sitioning accuracy is not good enough to spatially resolve 
any of these sub-components. However, the jet has a steeper 
spectrum than the core and its emission, which is already 
weak in the radio domain should be negligible at millimeter 
wavelengths. The relevant size of the background contin- 
uum, regarding millimeter absorption, should be < 4 mas 2 . 
Assuming a flat universe, with standard cosmological pa- 
rameters (H =70 km s _1 Mpc _1 , f2]yi=0.3, Oa=0.7), the an- 
gular size scale is 7.1 kpc/" at a redshift z — 0.68. The mm 
continuum emission of image A should therefore have an 
extent lower than ~ 30 pc 2 . This value, corrected from the 
filling factor, gives an estimate of the size of the molecular 
absorbing clouds. 

It is difficult to determine a direct value of the flux ra- 
tio Ja/Jb from our current data. Nevertheless, we have 
repeated the self-calibration procedure described in §2, by 
varying the flux ratio in the source model. The fitted posi- 
tions of the absorption, averaged over Vhel = —10 to +12 
km s , were then compared to the position of image A set 
by the source model (see Fig(3J) . Given the size of the contin- 
uum source, the average position of the absorption is consis- 
tent with the position of image A for Ja/Ib = 4.21^'q. We 
emphasize at this point that the average position derived in 
this way changes by less than ±20 mas when Ja/Ib varies 
from 3 to 6, i.e. by the range of possible Ja/Jb ratios, so 
that our conclusion that all the absorption arises in front of 
A is robust. Similarly, the self-calibration method that we 
have used depends little on the continuum source model, 
in particular on the distance between A and B, so that the 
uncertainties on this distance do not affect our results. 

Although indirect, our measurement of the flux ratio 
/a//b is the first at frequencies higher than 22 GHz. The 
value of 4.2 is slightly hig her than those mea s ured at 15 and 
22 GHz with the VLA (jO'Dea et all Il992t iPatnaik et ail 
ll993tlBTggs et al.lll999t ). and almost twice higher than those 
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Fig. 2. Relative intensity and position of the centroid of 
the HCO + (l-2) absorption as derived from the fit of a point 
like-source to the self-calibrated visibilities, for each veloc- 
ity channel, i) The intensity scale is zero at the continuum 
level and at —1 for full absorption. The profile of the ab- 
sorption, as fitted in FigJTJ is overlaid as a dashed line. 
ii, Hi ) Position offsets relative to the centroid of the con- 
tinuum, iv) Position of the four Gaussian velocity compo- 
nents, compared with the locations of the two lensed images 
of the quasar (grey disks, not to scale). 



observed around 2 GHz (|Mittal et al ll2006al). Th i s is con - 
sistent with the model developed bv lMittal et~ail (|2006bl ), 
where image A is obscured at radio frequencies by a HII 
region associated with the molecular cloud. We note that 
the Ja/Jb ratio might be affected by the time variability of 
the quasar, due to the time difference At between the tran- 
sit times of the light in A and B. The time delay, ho wever, 
is relatively short (At ~ 10 davs. lBiggs et al.lll999( ). mak- 
ing the chance of a factor of 1.5-2 variation of the quasar 
intensity in less than At rather unlikely. 

Adopting /a/Ib — 4.2 and assuming a filling factor f c 
= 1, we derive a maximum optical depth r max = 1.5 for 
component c, i.e. the peak opacity is large. Conversely, 
a lower limit to the filling factor f c can be derived by 
assuming r max — > 00 ; this yields f c > - J7. T he H 2 CO 
(IMenten fc Reidlll 996). NH 3 dHenkel et alJl2005l) and H 2 
(ICombes fc Wiklindlll997f) absorptions are probably caused 
by the same cloud as component c. For the other velocity 
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Table 1. Average positions of the different absorption components. For comparison, the positions of images A and B, 
as fixed in our continuum source model, are indicated. Offsets are given relatively to the barycenter of the continuum 
emission. 



Velocity 


Vhel t 


Absorption 


AV 1/2 T 


R.A. offset 


Dec. offset 


component 


(km s" 1 ) 


depth 1 * 


(km s" 1 ) 


(mas) 


(mas) 


a 


-6.3 


0.31 


4.7 


-84 ± 36 


+28 ± 75 


b 


-0.3 


0.32 


4.5 


-73 ± 24 


-28 ± 51 


c 


+5.9 


0.63 


4.5 


-60 ± 12 


-11 ± 26 


d 


+10.3 


0.30 


4.4 


-57 ± 38 


-27 ± 79 


a+b+c+d 


-10 to +12 






-61 ± 15 


-15 ± 32 


Image A 








-59 


-25 


Image B 








+250 


+103 



from the fit of the spectrum obtained in the compact configuration (cf. Fig[T]). 
with respect to the total continuum intensity. 
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Fig. 3. Average positions of the HCO + (l-2) absorption fea- 
tures obtained when varying the flux ratio /a//b (indicated 
above each position) in the source model used for the self- 
calibration. 



components, the filling factors and/or the peak opacities 
must be lower. A filling factor f c ~ 0.9 is consistent with the 
fact that image A, despite a high column density of absorb- 
ing gas, is strongly attenuated on optical V-band images, 
with respect to H-band images, but still visible, ft might 
also naturally explain why the optical distance AB (~ 317 
mas, lYork et al.ll2005h appears lower than in radio (~ 334 
mas) . The optical barycenter of image A should indeed be 
shifted closer to image B, if the absorption occurs mostly 
on the opposite border. 

By co- adding all our observations in the compact con- 
figuration of the PdBI at the frequency of the HCO + (l-2) 
line, we do not detect any new absorption feature, outside 
the Vhel = —10 to +12 km s _1 components just described, 
over the range Vhel = —200 to +1300 km s _1 , and this 
down to a level of 3cr = 3.7% of the continuum at a velocity 
resolution of 7.1 km s . As the intervening galaxy is seen 
face-on, it is unlikely that the difference in velocity between 
positions A and B, located at either sides of the center of 
the galaxy, is more than 200 km s -1 . HCO + is known to 
be easily detectable in a bsorption in the Gal axy, even in 
diffuse molecular clouds (|Lucas fc Lisztj fl996l) . Therefore, 
either the region intercepted by the light path associated 
with image B and located at ~ 400 pc from the center of 



the galaxy, is free or almost free of molecular gas (thco+ 
< 0.2), or the filling factor is low (/ c < 0.2). 

4. Conclusion 

Sensitive observations of the quasar B0218+357 in the 
HCO + (l-2) line, redshifted by z — 0.68, show at least four 
velocity components in absorption with velocities between 
Vhel = —6 to +10 km s _1 and widths ~ 5 km s -1 . By fit- 
ting the visibilities obtained from new observations with 
the very extended configuration of the Plateau de Burc 
Interferometer, we show unequivocally that all these com- 
ponents arise in front of the SW gravitational image (A) of 
the quasar. We see no other HCO + absorption components 
over a range Vhel = —200 to +1300 km s -1 , that could 
arise from image A, from the NE image B, or from the weak 
Einstein ring visible at radio frequencies. This implies cither 
a low average column density of HCO + (t < 0.2) in front 
of B and, presumably, a low average molecular hydrogen 
column density, or a low filling factor (/ c < 0.2). We derive 
a flux ratio Ja/Jb = 4.2+^q f° r the two main components, 
A and B, at 106 GHz, which is slightly higher than those 
derived between 15 and 22 GHz and almost twice larger 
than observed at frequencies around 2 GHz. This measure- 
ment, strictly speaking, applies only to Nov 8 & 9, 2005, 
the dates of our extended configuration observations. 
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